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ABSTRACT 

Turbulent, two-dimensional, hydrodynamic flows are characterized by the emergence of co- 
herent, long-lived vortices without a need to invoke special initial conditions. Vortices have 
the ability to sequester particles, with typical radii ~ 1 mm to ~ 10 cm, that are slightly 
decoupled from the gas. A generic feature of discs with surface density and effective temper- 
ature profiles that are decreasing, power-law functions of radial distance is that four vortex 
zones exist for a fixed particle size. In particular, two of the zones form an annulus at in- 
termediate radial distances within which small particles reside. Particle capture by vortices 
(~| , occurs on a dynamical time scale near and at the boundaries of this annulus. As the disc ages 

Oh' and the particles grow via coagulation, the size of the annulus shrinks. Older discs prefer to 

Q , capture smaller particles because the gas surface density decreases with time, a phenomenon 

i-{ • we term "vortex aging". More viscous, more dust-opaque and/or less massive discs can have 

^ , vortices that age faster and trap a broader range of particle sizes throughout the lifetime of 

the disc. Thus, how efficiently a disc retains its mass in solids depends on the relative time 

scales between coagulation and vortex aging. If vortices form in protoplanetary discs, they 

C^ ' are important in discs with typical masses and for particles that are likely to condense out of 

the protostellar nebula. Particle capture also occurs at distances relevant to planet formation. 

Future infrared, submillimetre and centimetre observations of grain opacity as a function of 

\^ ', radial distance will test the hypothesis that vortices serve as nurseries for particle growth in 

protoplanetary discs. 

'^ I Keywords: hydrodynamics -planets and satellites: formation 

o 

> 1 INTRODUCTION 



^\. ' The formation of ~ km-sized planetesimals from sub-micron-sized dust grains likely involves more than one physical process l lYoudirj 



*~i , I2OI O: Chiang & Youdin 2010). It is generally accept ed that particle growth in dusty, circumstellar discs is hierarchical, eventually forming 



planetesimals that are the building blocks of planets (lArmitagell2007Ll201(]|) . There is reasonable understanding of building small particles, 
but forming planetesimals from these particles remains mired in controversy. The initial stage of growth probably proceeds through the 
nucleation of sub-micron-sized dust grains from the primordial nebula, which then form the monomers of fractal dust aggregates up to 
^ 1 mm to ^ 10 cm sizes in > 10"^ yrs, beyond which growth is stalled by collisional bouncing and fragmentation JBlum & WurmluOOa : 
Zsom et alj|2010f) . In this regime, the particle dynamics and coagulation are described by Brownian motion and van der Waals forces. One 



of the best astrophysical pieces of evidence for grain growth to these sizes is the detection of 3.5 cm dust emission from the classical T Tauri 
star TW Hya (age ~ 5-10 Myr), located 56 pc away, which has a face-on circumstellar disc of radius 225 AU JWilner et alj|2005h . 

In standard models of protoplanetary discs, the gas pressure decreases radially outward. Gas in the disc then moves at sub-Keplerian 
speeds. Solid particles on Keplerian orbits experience a "head wind" with a velocity ~ 10^ cm s~ ^ — this head wind drags ~ 10 cm- to ~ 1 
m-sized particles into the central star on time scales of ~ 10-100 years ( IWeidenschillinall977ah . These time scales are much shorter than 



the characteristic time scale for these particles to collide and grow into larger particles that are unaffected by the head wind. lSafronovl ( Il969f) 
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Table 1. Fiducial parameter values 



Parameter(s) 


Description 


Adopted value(s) 


Mo 


initial disc mass 


0.01 Mq 


so 


initial disc outer radius 


20 AU 


M^,,T-,,R-, 


stellar mass, temperature, radius 


Mq,Tq,Rq 


Ps 


material/internal density 


3gcm-3 


KQ 


dust opacity 


1 cm2 g-l 


a 


viscosity parameter 


0.01 


7 


gas adiabatic index 


1.4 


M 


mean molecular weight of gas 


2.4 



and lGoldreich & Ward J1973[) suggested that this d ifficulty can be circumvented by dust settling into the midplane of the disc and triggering 
gravitational instability, but I Weidenschilling ([l980!) pointed out that turbulence generated by the settling impedes the process. 

Long-lived structures in the gas are a possible way of concentrating particles with sizes ~ 1 mm to ~ 10 cm and growing them 
to larger sizes. Such structures are usually high-pressure r egions in the gas that are capable of concentrating particles, which is a man- 
ifestation of Bernoulli's principlqj dKundu & Cohenl 120040 . Vortices are examples of long-lived structures — they are spiral, non-linear 
motions of fluid with closed streamlines. On Earth, ocean vortices have been observed to trap larval fish off the c oast of western Australi a 

WeizsackeJJl946l) . 



I Paterson et al.ll2008n . In astrophysical settings, the possible role of vortices in planet formation was suggested by 
based on the writings of Kant, in an article entitled Die Entstehung des Planetensyste ms ("The Origin of Planetary Systems"). Since then 
vortices have been suggested as possible nurseries for growth to > 1 m-sized particles f Ad ams & Watkins 1 19951; [Barge & Sommeria |l995 
TangaetalJl"9 96: Bracco et al. 1999; Godon & Livio 1999, 2000,; Johansen et al. 2004; Barr anco & Marcusl2005l: 



Fromang & Nelson 



Klahr & Bodenheimeill2006l ; llnaba & BargdbOOd ; IShen et alJl2006l ; lBodo et al.ll2007l ; lMamatsashvili & RicJ2009h . In stratified protoplan- 



etary discs, there are conceivable locations where the flow is turbulent and quasi-2D. An attractive feature of turbulent, 2D, hydrodynamic 
flows is that the fluid robustly s elf-organizes into large, coherent, long-lived vortices amidst a ba ckdrop of small eddies without a need to 
invoke special initial conditions JCarnevale et al.ll99ll ; IWeiss & McWilliamsll993l ; lTabelindl2002l) . Such a property arises from the fact that 
the so-called "vortex-stretching term" in the vorticity equation is absent in 2D, thereby allowing an inverse cascade of energy (and forward 
cascade of enstrophy). Turbulence thus becomes a friend and not a foe. If protoplanetary discs are capable of producing turbulent, quasi-2D 
flows, then these may seed large-scale vortices that may survive for many orbital time scales. Such optimis m should be tempered by the 



fact th at off-mi dplane vortices have not been observed in simulations of protoplanetary discs with dust settling ( IChiang|2008l ; I Johansen et al 
2009h 



The main question we are addressing in this study is: assuming vortices can be generated and sustained in discs, what sizes of particles 
do they capture, where are the capture locations and when does capture occurs? In we discuss/review the order-of-magnitude physics 
associated with protoplanetary discs, particle-gas interactions and vortices. In iJ3] we start with the simpliest case of a static, minimum mass 
solar nebula disc and show that there are generically four vortex zones within any disc with surface density and temperature profiles that are 
decreasing, power-law functions of r. In ij4l we consider the next level of sophistication, which is the case of an evolving, viscously-heated 
disc; we show that there are preferred locations and particle radii for vortex capture. In ij5] we generalize to the case of an evolving, viscous, 
irradiated disc. We demonstrate that discs which are able to both settle particles to their midplanes and capture them via vortices have upper 
limits to their masses that are consistent with most observed discs. We also show that the maximum particle radius for vortex capture in these 
discs is ~ 10 cm, independent of disc model and weakly dependent on stellar, disc and dust properties. In iJ6] we summarize our conclusions, 
discuss the open questions concerning the physics of vortices and describe the relevance of our results to observations. Table [T] lists the 
fiducial values adopted for the parameters of our models. 



2 ORDER-OF-MAGNITUDE PHYSICS: DISCS, PARTICLE-GAS INTERACTIONS AND VORTICES 

We consider vortices in an axisymmetric, viscous accretion disc with surface density E(r, t) and (effective) temperature T{r, t), where 
r is the radial distance from the star and t is the time. Disc material orbits with angular velocity fi{r) = (GM^/r'^Y'^ around a star with 
mass M-t,. The sound speed in a disc depends on T{r, t) and is typically 

"^ = V^^"""' ^ IiookJ ' ^'^ 

where 7 is the adiabatic index of the gas, fee is the Boltzmann constant, ttih is the mass of a hydrogen atom and we have ignored the radial 
dependence of temperature within the disc for simplicity. In the vertical direction, the gas is in hydrostatic equilibrium with vertical scale 



^ Bernoulli's piinciple states that in inviscid flows, a decrease in fluid velocity is accompanied by an increase in pressure. Hence, locations of maximum 
pressure are also locations of minimum velocity. 
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height 



Gas diffuses inward at a rate M 
where 



H 



£i 



0.03 AU 



T 



100 K 



1/2 



M0 



-1/2 



VAU/ 



(2) 



Zitv 



E, where ;/ is the viscosity. We adopt an "alpha-model" for the viscosity JShakura & Sunvaevll 19730 



ctCsH 



n 



(3) 



td ~n" 



: 0.2 yr 



/ r \3/2 
VAU/ 



There are three important time scales associated with the disc. The shortest is the dynamical (orbital) time scale, 

.Me, 

The disc cools on the local thermal timescale jPringldl 19811) . which is the ratio of the thermal energy of the gas (ECg) to the rate of viscous 
energy dissipation (j/Efi^). Using our expression for v (equation (3)), the cooling time is 

M^Y^^^ /^\3/2 

Jl^J VAU/ 

Thus, the disc cools on time scales much longer than the dynamical time. The viscous time scale is usually the longest characteristic time 
and measures the rate at which matter diffuses through the disc, 

^2 , ^ . 9. , ^ s 0. ''r/H'^ ^ ' "' ^ -1 / A/f \ "1/2 



icooi ~ a td ~ 20 yr 



VO.Ol/ 



(4) 



(5) 



tvis ~ — ~a (77) ^^^(77) icooi ~ 2 X lO"" yr 
It is apparent that ivis > icooi > td in typical discs JLvnden-Bell & Pringl3ll974l : |Pringlelll98lh . 



r/Hy (j^\-^ ( Ml.\ (JL.\ 
Too"; vooT/ Km^jJ vau/ 



3/2 



(6) 



In a Kep lerian disc, solving the radial momentum equation yields the velocity difference between a particle and its surrounding gas 
l Pringlelll98lh . For a thin disc, c. ^ uk where uk = (GA/*/?-)^". Defining inward radial drift to have a positive sign, 



At; 



2 4 4/ UK ' 



The quantities p and q are the power-law indices of the surface density and temperature profiles as functions of r, 

91nE ainT 



P = 



g = 



(7) 



(8) 



9 In r 9 In r 

From Equation l[7]l, a pressure gradient that declines radially outward produces an inward radial drift of particles. If p and q are sufficiently 
negative, then the particles do not drift. If vortices can provide this kind of environment, they will trap particles. Ignoring the dependence of 
Aw on p and q, we have 

— ^^~°-°ni5oKXuj [m;) ■ ^^^ 

Adachi et alj ( 119760 and IWeidenschillind ( Il977an derived 



For c. 



i;k ' VlOOKAU, 
1 km s^i (equation |T|), typical drift speeds are roughly 30 m s 



detailed expressions for the time scale of radial drift tdrift- As a rough estimate, tdrift ~ r/Au ~ 10^ yr (r/AU), longer than the dynamical 
or cooling timescales but much shorter than the viscous timescale or the expected dissipation time of the nebular gas (~ 10^-10^ yr). 
Computing tdrift more accur ately requires a careful consideration of the drag coefficients, which themselves depend on a d Adachi et alJ 
19761 : IWeidenschillin jl 1977ah . To facilitate discussion, we estimate an approximate radial drift time scale at r ~ 1 AU as 



tdrift ~ 10 yr 



(a/1 m)~i , a < 1 m, 
a/1 m, a > 1 m. 

In the absence of turbulence, the time scale for particles to settle to the midplane of the disc is JChiang & Goldreichll 19970 



500 yr 



M^ 

M^ 



-1/2 



/ r \3/2 

Iau/ 



(10) 



(11) 



apsQ """J \^io3gcm^2y ylmmSgcm^^^ 
where ps is the material/internal density of the particles and we have ignored the radial dependence of the surface density. Equation Jilt is 
effectively the minimum time for the particles to settle, a s vertical mixing may hold the particles aloft. Even turbulence as weak as a '-.^ 10"^ 
may frustrate the settling process (see lCuzzi et alj2008l and references therein). (If the particles are porous/fractal, the settling time may also 
increase.) Particles with radii a ^ 1 mm (tdrift > 10* yr) are well-coupled to the gas and settle somewhat slowly to the disc midplane. 
Particles with radii a ~ 1 mm to ^ 10 cm (tdrift ~ 10^-10* yr) are weakly-coupled and settle on time scales > 5-500 yr. The largest 
particles with radii 3> 1 m are uncoupled and do not "see" the gas. Thus, particles are considered "small" or "large" depending on the 
extent to which they are coupled to the gas, a concept we will develop later in the context of vortices. Small particles may grow up to some 
maximum particle size, typically < 10 cm, before starting to drift inward. We wish to show that vortices are capable of capturing particles in 
this size range as they drift radially inward through the disc. 

The typical growth time for small particles is shorter than the drift time. Particles larger than micron sizes are probably fractal aggregate 
structures (see Appendix IaTi. which simulations and laboratory experiments suggest take (Blum & Wurm 2008 : Zsom et al. 2010) 



Z 10 yr 



(12) 
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to for m. Once these particles form and start to settle out of the gas, the collisional time scale for an ensemble of solid particles is JHeng & Tremaine 
2O10h . 



^..xin-3v,-r^^ e^\( ^- y ( M^y"' ( j_\^'\ 



^-~S:»"'^^°" y'A lmm3gcm-3 jl,IFi7;^; (,A^J [W ' ^l^) 

where Es is the surface density of solids (typically 1% of E) and we have ignored gravitational focusing. If collisions result in the growth 
to larger particles, then particles with a ~ 1 mm to ~ 10 cm grow more rapidly (tc ~ 0.005-0.5 yr) than they drift through the gas 
(idrift ~ 10''~10^ yr). Larger particles with a ~ 1 m grow on timescales {^ 5 yr) comparable to their drift times. Thus, producing particles 
much larger than 1 m requires an environment to concentrate small particles in regions of larger local surface density where they can grow 
into much larger planetesimals which are safe from radial drift. Vortices provide this environment. 

Particles in the vici nity of a vortex experience centrifugal (due to vortex and not disc rotation) and Coriolis forces, directed outward 
and inward respectively JChavanijuOOOh . In polar coordinates centered on the vortex (_R, 0), the relevant terms in the equation for the radial 



acceleration (cPR/dt^) of the particle are: 

Centrifugal: R 



de^ ^ 



Coriohs: 2Q.R—. 
dt 



Only anti-cyclonic vortices (d9/dt < 0) direct particles toward the vortex centers; they also survive longer than cyclonic vortices JDavis et al. 



20001) . For net inward acceleration of the particle to exist, we must have \d6/dt\ < 2f2, a condition which is always fulfilled in practice. Two 



processes — viscous dissipation and orbital shear — limit the sizes of vortices. Viscous dissipation destroys vortices smaller than the viscous 
length scale, 

... ^ . 0.003 » (^^) (ii^) - (1^) -'' (i,)- , 



on time scales comparable to tvis (equation (6)), where iivor ~ 0.1 Cs jlnaba & BargduOOq) is the rotational speed of the vortex. Even 



in 



this case, vortices formed at this scale are many orders of magnitude larger than the particles conside red and can survive fo r many dynamical 
times. Keplerian shear inhibits the formation of circular structures larger than the shear length scale JGodon & LividI 19990 . 



^shca 



^^0.05AUr:^V7^r f^r^ (16) 



dr "'^^ \0.1cs) \Mq) KavJ 

Circular vortices forming at scales > Lshcar get elongated in the azimuthal direction, which allows them to survive longer. These estimates 
for Lvis and Lshcar show that vortices are large-scale phenomena with length scales L 3> a. 

In 3D, vortices are generally subjected to hydrodynamic instabilities that may destroy them on time scales shorter than ivis (see ij6]l. In 
2D, vortices may live for many dynamical times due to the absence of the vortex trapping term in the vorticity equati on. An inverse e nergy 
cascade exists in 2D turbulence, allowing large-scale vortices to emerge naturally from self-organization of the fluid (! Tabeling|2002i) . It is 



therefore much easier to make statements about vortices in 2D. The figure of merit for whether a fluid flow is 2D or 3D is the Froude number, 
which is a measure of how rapidly a fluid element responds, against the vortex flow, when vertically displaced in a con vectively stable fluid 



layer. At a given vertical height z in the disc, the fluid layer is effectively 2D if the Froude number is less than unity l lBarranco & Marcus 



20051) . It depends on the ratio of two quantities: the buoyancy time scale tev and the rotational period of the gas around the vortex tv 



Assuming vertical isothermality for simplicity, we get 



..v^l f I. 



As reasoned bv lBarranco & MarcusI 12003), the Keplerian shear is comparable in magnitude to the rotational rate, implying that ivor ~ 2-k/Q,. 



The Froude number is then 

^=^~-, (18) 

where we have ignored a factor of order unity related to the aspect ratio of the vortex. Near the disc mid-plane, J^ ^ 1 and the flow is 3D. 
A couple of vertical scale heights or more from the mid-plane, however, the flow already becomes quasi-2D — any turbulenc e generated in 



these fluid layers wil l inevitably seed vortices, but we again note that such off-midplane vortices are not seen in simulations dChianduOOa : 
Johansen et alj2009l) . 



ChavanisI 120001) considered the forces acting on a particle crossing a vortex. By seeking solutions of the form, 

X — Lcos (ujt) exp (— i/icap), 

(19) 
y = L sin (wt) exp (-t/tcap), 

where {x, y) are the Cartesian coordinates of the particle with the origin centered on the vortex, he was able to show that L' — L/x and 
cj — — 3f2/2(x — 1) — the particles follow ellipses of aspect ratio x ^ ^ and move with angular velocity oj. The eccentricity of the ellipse 
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is (x — 1) /x- The combined effect of the Coriohs force and drag make the particles drift towards the vortex center on a time scale: 

f 4x(x-l)^ (S 
, _ ) 3(x-2)(2x+l) \n 



4x(x-i)- r^) , ^>n 

2x(x-n /1^ ^ ^' (20) 



(x-3)(2x+ 



^(1)' ^<"' 



where ^ is the friction coefficient which we will describe shortly. The preceding expressions were derived by Chavanis (2000) under the 
assumption that the particl e motion within the vortex is deterministic, which becomes i nvalid if the vortex cores become strongly turbulent 



assumption ttiat trie particl e motion wittiin trie vortex is deterministic, wtiicti becomes i nvalid it ttie vortex cores become strongly turbulent 
I Lesur & PapaloizoJuOlOf) . Since x — 4 minimizes the capture time JChavanisluOOOl) , we will adopt this value for the rest of the paper 



whenever necessary and note that the numerical coefficients in equation J20l ( are both 8/3. The capture time depends on particle size and 
the properties of the gas through the friction coeffi cient f . Denoting the mean free path of the gas as A --^ 1 cm, particles are in the Epstein 
(Stokes) regime when a < 9A/4 (a > 9A/4) (e.g.. lAdachi et al.lll976l : IWeidenschillindll977ah . For "small" particles coupled closely to the 



ga s, the transition from the Epstein regime to the Stokes regime occurs at a disc radius r — r^ (see §3). Thus, we write the friction coefficient 



as 



Weidenschillinelll977al : lcuzzi et al.lll993l) 



t _ J J!^' Epstein regime (a < 9A/4 or r > re) , 



9EnA ^ Stokes regime (a > 9A/4 or r < r-c 



Solving for the vortex structure within the disc then requires specifying a disc model that in turn specifies E, T and A as functions of 
r and/or t, which we will explore in the subsequent sections (^ ij4]and ijsj- h is worth noting that Sl/^ is the commonly-used "Stokes 
number". If ^/Sl 2> 1, the particle is small and tightly coupled to the gas, implying that it is unabl e to be trapped a t the centre of the 
vortex. If ^/n <C 1 , the particle is large and does not "see" the vortex. Hence, as already noted by IChavanisI 1 120001) and simulated by 



Johansen et alj ( 120041) . vortices tend to "pick out" particles of a certain size, i.e., ^/fl ~ 1, such that capture occurs within an orbital period, 
tcap ~ Q~^. These particles are also the ones that attempt to settle towards the mid-plane of the disc after growing to large sizes via 
coagulation. This aerodynamic sorting of the particles in side vortices is consistent with evidence that meteorites in our Solar System are 



coagulation. Ihis aerodynamic sorting or the particles in side vortices is consistent with evidence that meteorites in our bolar System are 
composed of chondruleo of a similar size ( lHewinslll997 ). It is worthwhile to note that the expression for ^ in the Stokes regime only holds 



when the particle Reynolds number is less than unity. Equivalently, this requires 

- < stV (-) • ^22> 

"K 3 + q + 2p \a J 
Since Cs/vk <^ 1 for the discs we are considering, we expect the condition in equation J22b to be satisified. 

Finally, gravitational instability may be triggered if the concentrated mass density within the vortices exceeds the effective Roche density, 

where cc7 is a dimensionless factor accounting for the delay of gravitational collapse due to gas pressure. Setting zu — 1 gives the traditional 
mass density threshold for gravitational instability, which is about two orders of magnitude larger than the characteristic mass density of the 
gas within the nebula, po ~ 10~® g cm^"^. Traditionally, gravitational collapse occurs within 



where pc > pr is the mass density of the self-gravitating clump created by gravitational instability. Among others. ICuzzi et al.l ( 120081) point 
out that gas pressure and turbulence act to delay the onset of collapse, analogous to the role of ambipolar diffusion in star formation, estimate 
zu > 10^ and show instead that gravitational collapse occurs on a timescale, 

1 fpQCs^ ^2 _ ,nn.., f PO ^^ \ f P^ P^__2_\~^ n'i^ 

EG.eff ~ -r^ to ^ 400 yr — — ^ 3-— ^r 77—^? z^^ Z^^ ■ (25) 

Atv^ \ Psa J \10 ^ gem -^ 1 km s ^/ \ 10 ^ g cm •* 3 g cm ■^ 1 mm/ 

In summary, small particles (< 1 /im), strongly coupled to the gas, coagulate rapidly and attempt to settle to the midplane of the disc. 
Particles approaching ~ 1 mm to ~ 10 cm in size experience radial drag forces that direct them into the star. Vortex capture occurs on a 



— large-scale vortices therefore prov ide a way to trap these ] 
tices then enhances particle growth (Inaba & Barge 2006;. 



density within the vortices then enhances particle growth (Inaba & Barge 2006). Growing particles remain trapped within the vortices until 
they are large enough to decouple from the gas, when radial drag is negligible. Thus, vortices provide a natural way for small particles to 
grow into larger ones before they are dragged into the central star. 

This mechanism for trapping particles within the disc is an interesting alternative to streaming instabilities, the clumping of particles in a 
gaseous disc with pressure support (e.g., Goodman & Pindor 2000 ; Youdin & Goodman 2005 ; Johansen & Youdin 2007 ; Youdin & Johansenl 



20071 ; Ijohansen et al.ll2009l) . In streaming instabilities, efficient clumping requires dust-to-gas ratios of order unity near the disc midplane. 



Vortices require no special dust-to-gas ratio. Hence, vortices above the disc midplane and streaming instabilities in the midplane may provide 
complementary mechanisms for forming large planetesimals from small particles. 

To examine whether vortices can trap particles in protoplanetary discs, we consider several accretion disc models. To develop an initial 
picture of this process, we start with a standard, static protostellar disc where the surface density and temperature decrease radially outward. 

^ Chondrales are typically mm-sized particles found in chondrites, which are meteorites with near-solar compositions unaltered by heating processes. 
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To explore a broad range of discs, we consider masses and surface density gradients consistent witii observations of discs surrounding young 
stars. Our results suggest that four vortex zones generically exist in discs with surface density and temperature profiles that are decreasing, 
power-law functions of r. Two of these vortex zones form an annulus at intermediate distances from the star, within which particles are tightly 
coupled to the gas and are considered "small." Vortex capture of particles is optimal near and at the boundaries of this annulus. To explore 
the generality of these results, we then consider two models of evolving discs: a completely viscous disc (jQ and an irradiated, viscous disc 
(iU- Our analyses suggests that these basic conclusions are unaffected by the choice of disc model. 



3 VORTICES IN A STATIC, MMSN DISC 

Determining the locations where vortices trap part icles in gaseous discs involves specifying the surface dens ity and temperature profiles. 
In models of the minimum mass solar nebula (MMSN: IWeidenschillinslll977bl : lHavashilll98ll : ICuzzi et alj|l993l) . 



T - 



,'■0 



(26) 



At ro = 1 AU, the normalization values typically chosen are Eo — 1700 g cm~ and To ~ 280 K. In this section, we assume 7 = 1, 
i.e., the sound spe e d is i sothermal. Standard choices for the surface density and temperature indices are 1/2 ^ g ^ 3/4 and ^ p ^ 5/3 
I Cuzzi et al.lll993nMChav anis (2000) chose p — 3/2 and q — 1/2. Forp = 3/2, evaluating 2-k f 'Srdr from the approximate locations of 
Mercury to Neptune (0.22-35.5 AU) yields a total mass of about 0.01 Mq jWeidenschillind 1977a) . Repeating the same exercise with p = 
yields a mass of about 0.8 Mq . 

Submillimetre studies of young stellar objects show that the surface density drops off more slowly than predicted by the p = 3/2 
MMS N model, and that the distribution of disc masses is roughly log-normal with a typical mass ~ 0.01 Mq ( lAndrews & Williamsll2005 . 



20071) . Detection of continuum emissio n at 1.3 mm from 1 1 circumstellar discs around low- and intermediate-mass pre-main-sequence stars 
yields the constraint -0.8 < p < 0.8 Jlsella et al i ,2009): the disc mass ranges from 0.07 Mq (p = 0.8) to 9 Mq (p = -0.8) if we keep 



So = 1700 g cm~'^. In this ens emble, the most massive discs (~ 0.1 Mq) can produ ce giant planets by either core accretion jlda & LinI 
[2OO40 or gravitat ional instabili ty wafikov 2009, Kratter et al,| 20 id; Meru & BatduOlOl) . Lower mass discs may produce Pluto-mass objects 
and debris discs JHeng & Trem aine 2010; Kenvon & Bromlevl l2010l) . Inferred disc outer radii from partially-resolved submillimetre images 
of these ~ 1 Myr-old o bjects are typically ^ 200 A U; a viscous disc model with a ~ 0.01 matches the median spectral energy and surface 
brightness distributions JAndrews & Williamsll2007i ). Many discs with p < are dynamically unstable and may represent transient states 
that are less relevant for the vortex capture of particles, so we will focus on p Jj discs. 

Specifying the surface density and temperature profiles yields the vertical scale height and mean free path. 



H = Ho 



ro 



3/2-q/2 



A = 



2mw.H 
o-E 



= Ao 



''o 



3/2+p-g/2 



where 



Ao 



2mH / kjiTorl 



1/2 



(27) 



(28) 



,1/2 



For M^, = Mq, Ho ^ 0.05 



crEo VGA/^mn 

(7 ~ 2 X 10"^^ cm^ is the collision cross section of a hydrogen 
AU and Ao ~ 0.8 cm. The mass density is then po ~ Eo/iJo ~ 10^^ g cm^^. It follows that the pressure scales as P oc 7-~(p+9/2+3/2) 
Inward radial drift occurs when the pressure decreases outward, which requires (equation (T)) 

(g + 3) 



P> -- 



For a = 1/2, we need p > —7/4 for outward radial drift. 

ChavanisI J2000l) realized that there are four vortex zones within a gaseous disc (Figure [T](. Beyond a critical distance : 



(29) 



there 



is a transition from the Stokes (a > 9A/4) to the Epstein (a < 9A/4) regimes. The critical distance is bounded by two other transitional 
distances rin and rout. The annulus rin ^ r ^ rout defines a region within which particl es of a given size are tightly coupled to the gas. 
It is conceivable that coagulation takes place in such conditions ("e.g.. llnaba & Bargell2006r) . At the boundaries of the annulus (r = rin and 
r — rout), particles are captured by vortices within an orbital period (tcap ~ r2~^). Outside of the annulus (r < rin and r > rout), capture 
times become long and r adial drift domina tes. If particles grow rapidly within vortices, then r — r^ and r = rout are natural locations for 
planet formation. Indeed. Ichavanij ( I2OO0I) adopted a MMSN model and concluded that the Earth and Jupiter could have formed in-situ via 
core accretion at these locations. In a time-dependent disc, however, the locations and sizes of the vortex zones illustrated in Figure[T]evolve 
with time. Coagulation and vortex capture then need to work in tan dem to concentrate particles successfully. Before we embark on this task, 
we briefly review and generalize the static, MMSN results of lChavanis, C2000.) . 
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Figure 1. Schematic of different vortex zones witliin a gaseous disc. Small and large particles have ^/H ^ 1 and ^/f2 ■C 1, respectively (see text). 
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Figure 2. Transitional distances rin and Tout for a static, MMSN disc with q = 1/2 as a function of particle radius. The region with p ^ 0.8 is shaded, 
while the dotted and dot-dashed lines are for p = 3/2 and p = 1, respectively. The surface density and temperature normalization values are kept fixed at 
So = 1700 g cm-2 and Tq = 280 K, respectively. 



Upon specifying the disc model (equation 1261 '). we can compute the critical distance that separates the Epstein and Stokes regimes 



For a • 



\9XoJ 
Ao ~ 1 cm, Tc ~ 1 AU. In each regime, ^/f7 > 1 (small particles) when rin < r < r^ and r^ < r < Vo 

o X 2/(3-9) 



■ ro 



■ ro 



(l 

V9EoAo 



So 
2psa 



i/p 



(30) 



(31) 



Effectively, the disc can be divided into 4 sub-regions or zones as illustrated in Figure [T] Two of the zones form an annulus within which 
particles are considered small, even if the entire disc is populated with monodisperse particles (i.e., particles of the same size). Each zone has 
a different expression for the capture time: 



i 't^cap — ' 



27SoAo V'-o/ 



SSn ^' 



aAfl. 



3p 



(3-q)/2 

1 

(3-q)/2 



0. ( J-'] 



r < Tin, 

rin < r < Tc, 

Tc < r < Tout, 

r > rout- 



(32) 



In our idealized analysis, vortex capture of particles is optimal only at rin and rout • In reality, there will be a small range of distances around 
the transitional distances where vortex capture will still occur in about a dynamical time. We will see in ij4]that the boundaries of the annulus 
('"in ^ f ^ rout) evolve with time — specifically, the annulus shrinks as the disc ages. In ij5] the picture is complicated by the introduction 
of an additional transitional distance beyond which disc heating is dominated by stellar irradiation. 

The allowed ranges in the indices p and q imply a range of values in the surface density (Eo) and temperature (To) normalizations, but 
for clarity we keep the values previously described. We then plot rin and rout as functions of a in Figure |2] For clarity, we do not show re. 
The inner and outer transitional distances are independent of p and q, respectively. Since discs typically have a small range of values for q. 
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Tin for a given a is effectively constant for all discs; we keep q fixed at 1/2. As p decreases from 3/2, the outer transitional distance moves 
out for a fixed a. For each particle size, there is a set {rin, re, rout}. For all a values, there is a locus of sets for each combination of p and q. 
For illustration, we show rout for p = 3/2 and p = 1; we also shade the region for which p ^ 0.8. Decreasing the disc mass has no effect 
on rin, but causes rout to move inward. Thus, there is some value of a for which rin = rout. 

Particle capture via vortices is only possible in discs where rout > rin, which allows us to derive a maximum value for the particle 
radius, 



a < 



9 \ fcBTorgmH 



4cr / GM^, 



^] ^-(3+2p-,)/(3+4p-,)^ (33^ 



from which specializing to p = 3/2, q = 1/2 yields 

/ rr \ 3/17 / ,^ \ -3/17 X „ s 5/17 / s -11/17 

''<«"(5It<) (B im^) (ji^) ■ 

If less massive discs are considered, the largest particle that can be captured by vortices becomes smaller. For p = 0, rout is undefined. 
We demand that the particles settle to the midplane of the disc by a time 

t ^settle < imax, (35) 

where the dimensionless factor e ^ 1 accounts for an increase in the effective settling time due to vertical mixing. This constraint yields a 
minimum value for the particle radius, 

„ / 3 \ 1/2 / \ 3/2-p 

a > ' ' ' 



, (36) 



> 9 cm e 



\ 3/2-p / ,,r \ -1/2 / v^ \ / . 

r \ ' Ah\ ' En \ / p. t„ 



^roj \MqJ Vl700gcm-V V3gcm-3 lOyr, 

Equations l l33t and l l36t collectively yield a maximum value for Eo independent of a. In other words, we can estimate a maximum mass for a 
disc in which particles can both settle and be trapped by vortices. Let so and Sq denote the outer and inner disc radii, respectively; we further 
assume so ^ Sq . Combining equations ( I33t and ( 136b yields 

1/4 y„,, , (3+2p-9)/8p /, X (3+4p-q)/4p „so / ^ \ l-p-/3i 



Mo < 2'-^p+''-^^/*Ptv 



- — kBiomn 
Act 



(cm \'^''^^^~'*'I''^ ff \ WTip— y;/-ip rSQ / 



\ri ) \ . J - L U ' '-' '''' 



where /3i = (3 — 2p)(3 + 4p — <l)/8p. Specializing to p = 3/2 (/3i = 0) and q = 1/2 produces a more wieldy expression, 

/ e n \ 1/2 / T^ X 1/4 ^ ,, X 11/24 ^ \ 17/12 

-17/12 I So Ps \ / -to \ / JW* \ / tn 



" \20AU3gcm-3y VsSOKy/ \Mq J VlOyr,' 



This mass is comparable to the maximum disc mass observed in protostellar discs with ages of 1 Myr d Andrews & WilliamsluOOSl . 120071) . 
Note that there is no constraint on the disc mass when 

P=^. (39) 

For q = 1/2, this constraint occurs when p = 5/4. The mass constraint on a disc hosting particles that can both settle and be captured by 
vortices depends somewhat weakly on the disc properties and sensitively on the maximum time imposed for settling; f max may be interpreted 
as the radial drift time or even the gas dissipation time. Conversely, equation l l38b informs us that there should exist many discs in which the 
particles may be captured by vortices before having a chance to settle to the disc midplanes. In both scenarios, the settling/drifting particles 
may be sequestered by off-midplane vortices that are effectively 2D (see equation 1181 ). As demonstrated by the scaling dependence of e, the 
condition in equation ( 138b becomes more restrictive (i.e., lower maximum mass) when vertical mixing is present. 

Our static, MMSN model allows us to draw the following conclusions: in a disc where the surface density and temperature are power- 
law functions that decrease with radius, four vortex zones exist for a < 50 cm. Two of the vortex zones form an annulus at intermediate 
distances where particles of a given size are considered "small" even if the entire disc is populated with particles of the same size. Vortex 
capture occurs on a dynamical time scale near and at the boundaries of this annulus — capture occurs for particle sizes which coagulation 
likely produces in discs with masses comparable to those of observed discs. We next consider an evolving disc dominated by viscous heating, 
where the transitional distances discussed in this section and Figure[T]are allowed to evolve. 



4 VORTICES IN A VISCOUSLY-HEATED DISC 

ChambersI J2009t) considers an evolving disc with a structure determined solely by viscous heating. He solves equations for the vertical 



structure and energy balance to derive analytic expressions for the disc mass, radius, surface density, and temperature as functions of time. 
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The surface density and temperature are power-laws in radial distance and time, 

^'"^ (40) 

/ 8 \ 1/4 / ^ N -3/4 ^^"^ 

where so is now the initial value of the outer edge of the disc, kq is the opacity of the particles (assumed constant) and r = i/tvis- The 
coefficients Evis and Tyja are fun c tions of the initial disc mass and outer radius. We adopt sq = 20 AU as a reasonable starting point. We 
note that the models of IChambersI J2009l) do not con sider surface density enhancements at the snow line. 



The viscous time scale in the lChambersI 120091) model is 



167rci7fcBrvisEvis' 

where f^o = (GAf*/so)^", fi is the mean molecular weight of the gas and Mo is now the initial disc mass. Using the parameter values listed 
in Table[T] fvis ~ 4.4 x 10* yr. The normalization valuesijare 

7Mo 



lOTTsg ' 



/ 9 \ 1/3 

J- vis 



(42) 



\ MosB^iraYi 

where ctsb is the Stefan-Boltzmann constant. Using Table[T] Evis ~ 50 g cm~^ and Tvis ~ 27 K. 
From equation l l40t . the vertical disc height and gas mean free path are: 



r \ ,. . N_19/128 



9/8 

n = nA-) (i + r) 

n ZJ / \ 69/40 

2mH-H0 I r \ ,^ ^ _\361/640 



\^=^^^{- (1 + r 



(43) 



where 






The three transitional distances illustrated in Figure [T] are 

\9mHHoJ 

/2aaEvis 
V 9mnHo 

5/3 

rout = So 



/ /I 2 \ 8/9 

so(^^^) (l + r)i«/"\ 

/O V \ 40/69 

^c-^ofH^] (1 + .)-"/"-, (45) 



\2psaJ ^ ^ ' 



For convenience, we note that the temporal indices are approximately 0.13, -0.33 and -1.19. In an evolving disc, the positions of the four 
vortex zones change with time. The radius rin grows with time; the radii Tc and rout shrink with time. Thus, the region where particles of a 
given size are considered small (rin ^ r ^ rout) shrinks with time. Conversely, this region is larger for smaller particles at a given moment 
in time. 

Particles with ^/fi — 1 are preferentially captured by the vortices and by definition have Jlicap ~ 1 (equationl20t. These particles have 
radii of 

Again, for convenience we note that the temporal indices are approximately -0.074 and -0.71. 

As the disc ages, the vortices prefer to pick out smaller particles, a phenomenon we term "vortex aging". In Figure [3] we plot Oopt 
as a function of r at t = 10^, lO'^'^ and 10^ yrs, representi ng the possible range of coagulation time scales inferred from both theory and 
observations l IBlum & WurTnll2008l : IChiang & Youdinll2010f) . If our fiducial, evolving, viscous disc is capable of generating vortices, these 



vortices prefer capturing particles with radii ~ 1-10 cm (left panel of Figure[3}. If coagulation builds particles up to maximum sizes ~ 1 
mm, then at t = 10^ yrs vortex capture is only possible at r > 16 AU. 

Analogous to equation ( I33t in ^ the maximum particle radius imposed by the condition rout > rin is 

^ Equation (20) of lChamberj 420091) has a typographical error in the Jlp term. 
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Figure 3. Optimal radius of particles captured by vortices for evolving, viscous discs at t = 10'^, 10*'^ and 10*' yrs. The vertical lines are the respective snow 
lines (T = 170 K). Left panel: using the fiducial disc parameters (Table[T). Right panel: a disc that is more viscous, more dust-opaque and less massive than 
the fiducial case. Note that the vertical scales for the plots are different. 
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The minimum particle radius set by equation J35l > is 

9/10 



(47) 



(1 + ^) 
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By combining equations | |47I ( and ( I48t , it follows that discs which are able to accommodate both particle settling and vortex capture are 
subjected to a constraint on their masses, 



M < 0.28 Mq e 
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For convenience, we note that the temporal index is approximately 0.71. As in §3, the disc mass constraint depends weakly on the stellar and 
disc properties, but it is sensitive to the maximum time imposed for settling. For likely ranges of disc parameters, the maximum disc mass is 
consistent with the masses of observed discs. 

For vortices in discs to capture particles, they need to age on a time scale comparable to that required for coagulation, and to attain 
capture radii comparable to the maximum particle size attainable via coagulation. Lower values of Oopt are obtained when discs are generally 
more viscous, more dust-opaque or less massive. We have examined each of these cases separately, but do not show them; instead, we plot 
the extreme example of a more viscous, more dust-opaque and less massive disc in the right panel of Figure[3] Note that the vertical scales 
in the left and right panels are different. In this case, it is clear that Oopt ~ 1 mm to ~ 10 cm for t ~ lO'^-lO^ yrs. Almost the entire disc 
is in the Epstein regime, where the vortices age more rapidly. Our results may be invalid when r <C 1 AU because non-linear Stokes drag 
becomes important, which occurs when the gas Reynolds number, 

2aaEv 



7^' 
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Sn r 
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mn 
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greatly exceeds unity. It follows that non-linear Stokes drag sets in when 

V \ 5/8 
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5 VORTICES IN A VISCOUS, IRRADIATED DISC 

ChambersI ( l2009h also considers the case of an evolving, viscous, irradiated disc. Whether a disc has an irradiated outer region depends 
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upon comparing two characteristic temperatures, 
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lOTTsg 
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(52) 



Jj \Tc/ \S0 

where Tc = GAI^^niH/kBR*, T* is the stellar effective temperature and _R* is the stellar radius. Using Table[T] T^is ~ 27 K and Trad ~ 38 
K; since Trad > Tvis, the fiducial outer disc is dominated by stellar irradiation rather than viscous heating. However, if Afo, no or a are 
increased by an order of magnitude, we get Trad < Tvis — the entire disc is initially dominated by viscous heating. For clarity, we first 
consider the case where the disc has an initial irradiated region; we subsequently consider a disc where the irradiated region develops at a 
later time. 



5.1 With Initial Irradiated Region (Trad > Tvis) 
The characteristic surface densities are 



2_/vis — ^r 



4/5 



^rad — 
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287rs2 
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~98 \T^ 



52/33' 



(53) 



Using Table[T] Evis « 54 g cm~^ and Er, 
respectively; our fiducial disc now has ivi. 
The surface density profile is 



Requiring the surface density to be continuous yields the transitional radius between the viscous and irradiated regions, 



=:i 41 g cm ^. When Trad > Tvia, we replace Evia and Tvis in equation flTt by Erad and Trad, 
3.8 X 10" yrs. 



(54) 




n 



So 



70/33 



(1 + r 



-133/132 



(55) 



For our fiducial disc, the initial value of rt is about 1 1 AU. 

This formulation leads to minor discontinuities in the radial temperature. In the viscous regime, the disc is not vertically isothermal. 
Thus, the effective temperature is lower than the midplane temperature. In the irradiated regime, the assumption of vertical isothermality 
requires a discontinuity in one temperature at the boundary between the two regimes. .Chambers QQQS) adopts a model where the midplane 
temperature, 

Trrid=Trad<^ ^'\ ^3°// (56) 

ife) ' 

is continuous by construction, but the effective temperature, 

Ife) ' 

is then discontinuous at r = rt. Despite this lack of continuity, the jump in the effective temperature at the boundary is fairly small. Thus, 
our eventual estimates for aopt near r = rt should be correct to within a factor of a few. 
Knowledge of the effective temperature allows us to derive the vertical scale height. 
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(58) 



(59) 
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Figure 4. Transitional distances for an evolving, viscous, in'adiated disc for a = 1 cm particles and adopting tiducial parameters (Table[TJ. The yellow shaded 
region is the region within which a = 1 cm particles are considered small. In the pink shaded region, the small particles are in the region of the disc where 
heating is dominated by stellar irradiation rather than viscosity. 



It follows that the gas mean free path is 
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The critical transitional distance between the Stokes and Epstein regions now has two expressions, 
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(viscous), 
(irradiated). 



(60) 



(61) 



For the range of parameters considered in this study, we always have re < rt at t — 0, implying that the discs begin with Vc in the viscous 
regime. Since rt decreases more strongly with time (eauationl55b. it is possible for some discs with a monodisperse population of particles to 
have Tc in the irradiated regime at a later stage of evolution. For our fiducial disc, rc(f = 0) ~ 0.3 AU. The other transitional distances are 



and 



r\n = So 



fout = So 



^9mH-Ho 
/ 4psa^cr \ 



8/9 



(1 + r 



19/144 



7/9 



(, 2p"'a ) ( 



1 + r) 



-19/16 



2ps fi 



14/15 



(1+^) 



-133/120 



(viscous), 
(irradiated), 

(viscous), 
(irradiated). 



(62) 



(63) 



Although we have shown both cases for completeness, we typically have nn(f — 0) and rout(i = 0) in the viscous and irradiated regimes, 
respectively; their fiducial values are 0.005 AU and 120 AU. 

Formally, there are four different scaling relations for the optimal capture radius, depending on whether rt > Vc or rt < re: 

'^/^^fe)'''^! + -)"'''''' ^<-c<rtorr<rt<re, 

V 4,psa \soJ 

A generalization of the previous expressions for Oopt (equation|46t is that there are now three scaling relations for the optimal capture radius 
for any given disc. In essence, the first and fourth equations in l |64| l are always in use; whether the second or third equation is used depends 
on if re is in the viscous or irradiated regime. 



Oopt 



re <r < rt, 
rt < r < re, 
re < rt < r or rt < re < r. 



(64) 
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Figure 5. Optimal radius of particles captured by vortices for evolving, viscous, irradiated discs at t = 10"^, lO'*'^ and lO'^ yrs. The vertical lines are the 
respective snow lines (T = 170 K). Left panel: using the fiducial disc parameters (Table[TJ. Right panel: a disc that is more viscous, more dust-opaque and 
less massive than the fiducial case. Note that the vertical scales for the plots are different. The kinks in the curves are artefacts of the discontinuous nature of 
the effective temperature in the viscous, in'adiated disc model (see text). 



In Figure|4l we show the evolution of the four transitional distances for our fiducial disc and a = 1 cm particles. We have rout > sq 
initially and even after one viscous time scale, i.e., the disc has no Epstein region for large particles. As the disc evolves, the region within 
which a = 1 cm particles are considered small (i.e., they are tightly coupled to the gas) shrinks, although it does not disappear before the gas 
dissipates (r ~ 100-1000). (The region bounded by rin ^ r ^ rout is generally larger for smaller values of a.) The outer disc (r > 11 AU) 
is initially in the irradiated regime — as the disc ages, stellar irradiation becomes increasingly important as expected. At r ~ 200, the disc 
transitions from re < rt to re > rt, consistent with our earlier statement about monodisperse discs evolving to having re in the irradiated 
regime. Monodisperse discs with lower values of Mq or a require less viscous times to evolve to the re > rt stage. 



5.2 Irradiated Region Develops Later (Trad < Tvis) 

Discs with Trad < Tvis begin their lives entirely dominated by viscous heating as described in ^ However, at a time 
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stellar irradiation in the outer disc becomes important. Instead of Mo and so, the irradiated region is initiated with 



Ml = Mo 



-^ rad 
-^ vis 



21/73 



Si = So 



-^ vis 
-^ rad 



42/73 



(65) 



(66) 



Upon estimating rappoar. Mi and si, one then revises the values for the characteristic temperatures, 

'27«:oa7itB!^o\^^^ f 7A/o V'^ f 27Koa'ykBQiY^^ f IK'h 



J- vis 



64asBM"T-H 



lOTTsg 



Trad = r* 



4x1/4 .y. 1/7 .^^ 



3/7 



T, 



64(TSBM"^H 



4\ 1/4 /y \ 1/7 



2/3 



Te 



.7 1 V Tc y V So y V 7 

where fJi = (GM*/sf )'^' ■^. Other quantities are similarly revised or defined: 
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(68) 



The expressions for H and A (equation 1431 ) are now replaced by those in equations {58} and l l60t but with all of the revised quantities 
substituted. The quantities rin, re, rout, rt and Oopt are also modified/introduced in a similar spirit. 
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5.3 Optimal Capture Radii 

Following !J4]and Figure [5] we plot Oopt for evolving, viscous, irradiated discs in Figure |5] As before, we show in the left panel the 
fiducial case. The kinks in the curves are artefacts of the effective temperature profile being discontinuous at r = rt, as previously noted in 
i]5.1l Nevertheless, it allows us to conclude that 1 mm < flopt ;$ 10 cm, similar to our conclusions in ij4] In the right panel of Figure |5] we 
again demonstrate that more viscous, more dust-opaque and less massive discs are more amenable to sequestering particles in vortices. Both 
the discs shown in the left and right panels of Figure[5]have an initial irradiated region (Trad > T^is)- 

While the exact value of a^pt at a given value of r differs slightly between the viscous (@ and viscous, irradiated models, our general 
conclusions are the same. 

5.4 Particle Size and Disc Mass Constraints 

We generalize equations {47}, J48t and ( |49t to obtain constraints on viscous, irradiated discs that are capable of having particles settle 
to the disc midplane and also capturing them via vortices. The minimum particle radii obtained from demanding e igcttic < imax are 



ps^Qt 



^^^^^1^) (l + r)-«/^«, r>n, 



which can be rewritten as 



-"Z**" r<n, 



{mJ l20AU; Ugcm-aiOyrj | (-f^^ (,,|-j^) (1 + .)-"/- , r > n. 

For the algebra to be tractable, we make the approximation that Evis ~ TMo/Wttsq only in the expression for T^is (equation 1521 ). While 
this is formally not self-consistent (see equation 1531 ), the correction made to Tvis is small. It then follows that the maximum particle radius 
derived from imposing the constraint rout > '"in is 

/^ \ -376/3255 ^^ .352/1085/ J]' \ 242/465 / j, x 32/1085 / ^ x 16/1085 

"<'''^"Im;J (soau) Isoi^j [%) UJ (,i) 

/ fl N -4/105 / J X 76/465 . a N 16/465 f p, \ "^^''^^ f Kp \ '^^'' 95/248 

""{2.4) \1a) lo.OlJ [Sgcm-^J \lcm^g-^J ^^^^ 

in the viscous region of the disc (r < rt). While the maximum radius is the same as before (equation 1471 ). the scaling dependences are 
different because stellar irradiation modifies the temperatures and scale height. In the irradiated region of the disc (r > rt), the maximum 
particle radius is 

/ M^ \ "^/^^ / so \ 213/560 / Evis \ ''^^'"' fn\ "'^''"" / R ^ "'^'''*° 

a < 41 cm 



MqJ V20AUy V50gcm-2y yy^y y^^^^ 

71120 .^. 1/10 .^. 1/10 f Ps \ ""^'' f ^0 \ '^"' . .-57/128 

.2.47 V1.4/ VO.Ol/ V3gcm-3y \l cm'^ g-^ J ^ ' 

We next derive separate constraints on the disc mass in the viscous (A/vis) and irradiated (Afrad) regions of the disc. Combining 
equations {69} and J71t . we get 

/„ /„/ \ 1063/2230 / »^ \ fe , „ .-1409/3122 /T \ 96/1561 /„ s 48/1561 



V20AU/ 



100 / \MqI V20AU/ \TqI \Rq, 

._^. -124/1561 .^^.16/223 / p, V^''^^^^ ( np V^''' f t^V''^''' 2451/3568 

V2.4/ Vl.40.01/ VSgcm-Sy/ I^lcm2g-V UOyry' ^ ^^ 

where /32 ~ 186056/232069 ~ 0.80 and s'p <^ sp again denotes the inner disc radius. For the irradiated region, combining equations (|69} 
and l l72t yields 

-2/7 , „„ ....7/14 /T^ \2/7 / r, \ 1/7 



0.15 \M^) V20AU/ \%j) \R^ 



M.ad < 0.025 Me .-/^^(i;^] (^) 



V2.4/ V3ecm-3y VlOyry ^ ^ ' 



(74) 



19/16 



^3gcm~3y Vioy/ 

where the function ijf) is defined as 

<^(rt/so) = l-(rt/so)"'"°. (75) 

For n/sp = 0.25-0.75, « 0.29-0.07; (/)(0.5) « 0.15. We note that 0.025 Mq ^ 26Mj, where Mj is the mass of Jupiter. The main 
differences from the viscous region are that the dependence on so is somewhat stronger and there are no dependences on a or Kp. As 
expected, the disc mass in the irradiated region makes only a modest (~ 5%) contribution to the overall disc mass M = Mvia + Afrad. In all 
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of the models (static MMSN, viscous and viscous, irradiated discs), the maximum disc mass has strong dependences on imax and e with the 
power-law indices spanning about ±1.4-2.2; the dependences on stellar, disc and dust properties are somewhat weak. 

We conclude that the maximum radius for particle trapping via vortices is ~ 10 cm, independent of disc model and with somewhat weak 
dependences on stellar, disc and dust properties. We also conclude that a disc that is able to both settle particles to its midplane (in 10 yr) 
and capture them via vortices has a mass that is at most ~ 0.1 Mq, consistent with the masses of most observed discs. This mass threshold 
decreases if vertical mixing is present in the discs to hold the particles aloft for longer than fsettie. Nevertheless, particles may drift and/or 
settle from off-midplane locations and be captured by vortices in about a dynamical time. 



6 DISCUSSION 

6.1 Summary 

We have examined particle trapping by vortices in evolving, viscous and/or irradiated discs as a function of radial distance, initial disc 
mass, dust opacity and viscosity. The salient points of our study are: 

(i) If the surface density and effective temperature of protoplanetary discs are decreasing, power-law functions of the distance from the 
star, then all discs contain four vortex zones for a fixed particle size (Figures[T]and|2}. 

(ii) There is an annulus at intermediate distances from the star where small particles reside and may grow via coagulation. The size of this 
annulus decreases with time and also shrinks as the particles grow (Figure 0. Vortex capture is optimal near and at the boundaries of this 
annulus, and occurs within an orbital period. 

(iii) The optimal capture radii (~ 1 mm to ~ 10 cm; Figures[3]and[5]l of particles are comparable to or smaller than the sizes of particles 
(~ 1 m) that drift fastest through the protoplanetary disc. The enhanced particle concentration within the vortices may make the characteristic 
time for particle growth shorter than that for radial drift. 

(iv) Vortices in older discs prefer to capture smaller particles, a phenomenon we term "vortex aging" (Figures |3] and |5]l. If coagulation 
between sub-micron particles mixed with the gas can only produce small (< 1 mm) particles, then vortices can capture them throughout the 
lifetime of the disc. However, if coagulation manufactures larger (~ 10 cm) particles when the disc is young, vortices must form early to 
capture them. This metaphorical dance between coagulation and vortex aging determines how efficiently vortices help the disc to retain its 
mass in solids. 

(v) More viscous, more dust-opaque and/or less massive discs can have vortices that trap a broader range of particle sizes throughout the 
lifetime of the disc. While the coagulation of grains with sizes < 1 mm needs to be more in synch with the evolution of more massive discs, 
such discs are also expected to grow grains to larger sizes more rapidly. 

(vi) The maximum size of particle that can be trapped by vortices is ~ 10 cm, independent of disc model and weakly dependent on stellar, 
disc and dust properties. Discs where particles settle to the midplane (within > 10 yr) and are sequestered in vortices have upper limits 
to their masses (> 0.1 A/q) that are consistent with those of most observed discs. If vertical mixing is present (e.g., via turbulence), the 
maximum masses can be much smaller. 



6.2 The Physics of Vortices: Open Questions 

Many open questions remain concerning the microphysics of vortices. While we have shown that particles can be gathered by vortices, 
the outcome of these captures is uncertain. Vortices may concent rate enough ma ss to enhan ce collisions rates b y an order of magnitud e 



(or more) or to trigger local gravitational instabilities JAdams & W atkins 1995; Godon & Liv io 1999, 2000; Kla hr & Bodenheimenl2006l) 



Both paths leads to the formatio n of self-gravitating objects which will not drift through the disc. However, at least some of the published 
simulations are run in 2D (e.g.. lOavis et alj|200(l ; llnaba & Bargdl2006l ; iLvra et alj|2009l) — in the absence of viscosity or particles, such 



vortices live forever, implying that the mass of the planetesimal, embyro or planet formed depends either on the time the simulation is 
executed or the breakup of the vortex by the non-linear feedback of the concentrated particles. It is more likely that if the centres of vortices 
are relatively quiescent, they then serve as n urseries for c oagulation to occur between somewhat larger particles. 

Vortex formation is also uncertain. Les ur & Papalo izou (2010) show that the subcritical baroclinic instability (SBI) is a plausible way 
of seeding vortices. Discs develop this non-linear instability when they are (radially) convectively unstable, have non-negligible thermal 
diffusion, and are subjected to finite vorticity perturbations (~ 0.1). The associated Reynolds number for the shearing box simulations in 
this study is 7?. ~ 10'' ; because the threshold vorticity amplitude for invoking the SBI decreases with increasing TZ, they speculate that the 
threshold amplitude could be very small (and possibly sub-sonic) in realistic discs. 

Another possibili ty for generating vort ices is via the (linear) Rossby wave instability ^Lovelace et al.*4999;^ Vamiere & Taggenl2006[) , 



which was invoked bv llnaba & Bargd yOOg) to consider 2D vortices in protoplanetary discs. Inaba & Barge ( 2006) found that the formation 



of vortices via the Rossby wave instability critically depends on the amplitude and width of an initial density bump placed within the disc. This 
bump appears to be most unstable to perturbations with an azimuthal mode number of 5 (i.e., "m, = 5" perturbations). iMeheut et al.l J2O10l) 



performed 3D simulations of stratified discs and concluded that strong and persistent vortices emerge out of the flow via the Rossby wave 
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instability. IPavis et alj ( 1200(1) have noted that Rossby waves are only supported by flows with non-vanishing vorticity gradients, implying 
that they are relevant mostly in incompressible, inviscid flows. ^^^^^_^^^^^^ ^_^ 

The subject of vortex survivability is mired in deeper controversy. iLesur & Papaloizoul ( 120091) subject vortices embedded in a shearing 



sheet to 3D perturbations; only vortices with 4 < g < 6 survive the elliptical instability in unstratified discs. This stability region van- 
ishes when st ratified discs are considered. By contrast, Lithwickl ( 120091 ) asserts that weak vortices (g 2> 1) can survive in quasi-2D flows. 



Lesur & Papa loizou ( 20 1 0) find that vortices develop bursts of turbulence in their cores before surviving as weaker vortices ; the SBI amplifies 



the vortices and the cycle restarts. If turbulent vortex cores are generic and ubiquitous phenomena in realistic discs, then they may pose a 
setback to using vortices as mechanisms for concentrating particles. 

Many collective properties of vortices in turbulent, 2D, hydrodynamic flows remain poorly understood. Among these is the "universal 
decay theory", which is the empirical observation that the vortex density, radius, velocity, mean separation between vortices, enstrophy and 
kurtosis can be approximated by power laws of time parametrized by a single parameter — a first-principles explanation is still being sought 
([labeling 2002). Statistical theories predict vortices to always ultimately merge, in contrast to experiments that show that above a critical 
separation, a pair of vortices may remain separated for many dynamical times (lTabelin3l2002h . ^^^^^—^^^^^^^.^ ^^^^^_^^^^^_ 

Despite these uncertainties, vortices provide an interesting alternative to streaming instabilities (e.g.. lGoodman & Pindoii2000l : lYoudin & GoodmanI 



12005: Johansen & Youdin 2007; Youdin & Johansen 2007; Johansen et al. 2009), which require dust-to-gas ratios to approach of order unity 
(presumably near the disc midplane). Vortices require no special dust-to-gas ratio to trap particles and may exist off the disc midplane. In 
both cases, particles with ^/f2 ~ 1 are captured and the capture size decreases as the gas in the disc dissipates. Therefore, vortex trapping 
and streaming instabilities may provide complementary mechanisms for particle concentration and/or growth. A key question to study and 
explore is the size distribution of planetesimals produced by each mechanism, since this might have an impact on the types of (exo)planets 
produced in a system. Pr eliminary analyses of the i nitial size distributio n of Solar System plan etesimals suggest a broad range of values, 
ranging from ~ 1-10 km JKenvon & Bromlevlboid) to ~ 100-1000 km jMorbidelli etai]|2009l) . 



6.3 Observational Relevance 

Our study makes falsifiable predictions about the size of particles concentrated as a function of distance from the star. Particles smaller 
or larger than the optimal radii for vortex capture will be uniformly distributed throughout the disc, while those with sizes ~ 1 mm to ~ 10 
cm will be concentrated near their respective transitional distances rin and rout (Figure [T). 

If particle concentration leads to the growth of larger particles, these particles will then decouple from the gas and migrate out of the 
vortices. The decoupling will be stronger in the inner regions of the disc as ^/Q. oc a~^ (instead o f oc a~^ in the outer regions). In our Solar 



System, the smallest constitu ents (chondmles) o f ~ 100 km-sized asteroids have sizes ~ 1 mm jHewinsll 19971) . comparable to the sizes of 



particles trapped by vortices. ICuzzi et alj ( |2008» (and references therein) have pointed out that there is a spread of about I Myr between the 



formation times of the oldest and youngest objects in the same meteorite, implying that particle growth was fairly inefficient. This inference 
is in turn consistent with the limited temporal windows for particle trapping implied by vortex aging. Another interesting property of ~ 10- 
100 km-sized asteroids is that many of them are formed from a physically and chemically homogeneous mix of particles of a similar size, 
consistent with the aerodynamic sorting property of vortices. 

In extrasolar settings, observations characterizing grain opacity (and hence grain growth) in protoplanetary discs can quantify particle 
populations as fu nctions of radial d istance, but these are still nascent and have only been accomplished for a small number of objects, 
e.g., HD 163296 ( iNatta et al.ll2007l ). Therefore, the hypothesis that vortices serve as nurseries for particle growth will need to be tested 



by future infrared, submillimetre and centimetre observations of protoplanetary discs, and may lead to constraints on the time scales for 
grain coagulation (Figures[3]and|5}. For example, the Atacama Large Millimeter Array (ALMA), which operates at wavelengths between 0.3 
and 3.6 mm, might be able to detect disc features associated with vortex capture. Coagulation-fragmentation simulations can subsequently 
be performed to co nvert the observed fluxes and a values into constraints on the grain properties (e.g., mass/size distribution, porosity; 
Birnstiel et aklboioh . 



Even if vortices do not survive long enough to create self-gravitating structures, they will certainly play a strong role in redistributing 
matter throughout the disc. Whether the redistribution of matter by the vortices plays any significant role in the eventual formation of 
planetesimals is unknown. An implication of a size-dependent redistribution of matter is that if the amount of electric charge carried by a 
particle is proportional to its size, then charge separation in protoplanetary discs may be a fairly common phenomenon. 

Our study has shown that if vortices form in protoplanetary discs, they are important in discs with typical masses and for particles that 
are likely to condense out of the protostellar nebula. The capture of particles also occurs at distances relevant to planet formation. With this 
study, we hope to (re)ignite the debate in connecting the microphysics of vortices with the global properties of protoplanetary discs, as the 
first step towards understanding the efficiency of planetesimal — and eventually planet — formation. 
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APPENDIX A: COAGULATION OF MONOMERS IN A VISCOUS, IRRADIATED DISC 

One can estimate if binary collisions between two spherical dust grains will result in coagulation. When two elastic spheres collid e with 
relative velocity Wcoi, the characteristic binding energy associated with the compressed surfaces of contact is dChianp & Youdinll2010r) 



i^Mnd ~ 5a' (^y^', (Al) 

where 5 is the surface tension from unsaturated bonds and tp is Young's modulus, a measure of the "stiffness" of elastic material. Requiring 
the collisional energy to be less than the binding energy yields an upper limit for the relative velocity, 

/ 3(5 \ _i/2 , -1/3 , . -, 

Hence, two colliding particles will stick if Au < Ustick, where the velocity difference between a particle and its surrounding gas is given by 
equation Q. The maximum radii of these particles are 

•"•(&)""(^)""°(^)""°e + ^)"'" -- ,„, 

6/5 / \ -3/35 \^^) 

^397fcB7 \Ja) r > n. 

It is important to note that equation l lA3b makes no statement about how long it takes for coagulation to occ ur. ^_^ 

The values to adopt for Young's modulus and the surface tension warrant some discussion. Table 3 of IChokshi et al.l (19931) states that 




Tp — 7x10 erg cm~ for ice; they also list Young's modulus as 10 erg cm~ for graphite and 2 x 10 erg cm~' for iron. Terrestrial rock 
typically has ip ~ lO'^^-lO^^ erg cm~^. Given the uncertainties associated with dust composition and chemistry, we consider a generous 
range of ip ^ 10^*^-10^^ erg cm~^. However, since amax oc ip~^'^ , a variation of two orders of magnitud e in ip corresponds to a change 



of only a factor of about 6 in flmax. A somewhat larger range of uncertainty exists for the surface tension. IChokshi et alj ( 119931) estimate 
5 = 75, 370 and 3000 erg cm~^ for graphite, ice and iron, respectively — we then consider 5 ~ 10-10^ erg cm~^, which corresponds to 
two orders of magnitude of uncertainty in amax. Collectively, we get 1.6 x 10~* < 5ip^'^'^ < 0.10. With our fiducial disc, we find that 
maxjamax} ~ 0.01-0.1 /im at r ~ 5-6 AU for t < 1 Myr. Monomers of such sizes a ggregate to forrn fractal dust grains, whi ch are able 
grow to sizes ~ 1 mm because they possess internal modes of kinetic energy dissipation l IBlum & Wurmll2008l : IZsom et al.ll2010l) . 
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